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Abstract 

We explain the theory behind our fiber-optical analogue of the event hori- 
zon and present the experiment in detail. 
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A Theory 



In this appendix we describe the theory behind our fiber-optical analogue of the 
event horizon p]. After a brief summary of Nonlinear Fiber Optics we show how 
an optical pulse establishes a moving medium and how this medium corresponds to 
a space-time geometry. This point of view is relatively unusual in the fiber optics 
community, but, as we demonstrate, it is completely consistent with the established 
knowledge of this field |2J. We then take advantage of our approach in quantizing the 
electromagnetic field of a probe in the presence of a pulse. We describe the classical 
and quantum physics of fiber-optical horizons, explaining the classical frequency 
shifting and the quantum Hawking effect with as few assumptions on the physics 
and the prior knowledge of the reader as possible without going into excessive detail. 

A.l Kerr nonlinearity 

Optical fibers [21 E] are metamaterials — materials with optical properties that are 
dominated by their structure — for the following reason: an optical fiber confines 
light to narrow transversal regions, usually at the core of the fiber. The electromag- 
netic field in the transversal plane depends on the frequency and polarization of the 
light in relation to the transversal structure of the fiber. In this way [2], the waveg- 
uide gives rise to an effective dispersion [1] and birefringence [1] that supersedes the 
natural optical dispersion of silica glass. The confinement of the light in the fiber 
also enhances the intensity over the fiber length such that the natural optical non- 
linearity of glass becomes relevant, in particular the Kerr and Raman nonlinearity 
caused, respectively, by the electronic and molecular optical response of glass [2]. 
For making artificial event horizons, one can exploit the custom-designed dispersion 
and birefringence and the enhanced nonlinearity of optical fibers, in particular of 
microstructured fibers (also called photonic-crystal fibers) [3]. 

Let us first describe the optical properties of glass. In the laboratory frame, 
r = (x, y, z) describes the position in Cartesian coordinates and t the time. We 
use the operator V to denote spatial derivatives and dt for time derivatives. We 
characterize the intrinsic linear optical response of the glass fiber by the susceptibility 
profile Xg(^;r) and denote the nonlinear polarization vector of the medium by P. 
The vector E of the electric field strength obeys the wave equation in SI units. 



with ^0 being the electric permeability of the vacuum. We assume a lossless medium 
where Xg is real. In isotropic materials such as glass, no second-order nonlinearity 
exists on symmetry grounds [6]. The lowest-order nonlinearity is proportional to 
the response from the products of three electric- field components. Such an effect, 
called Kerr nonlinearity [2] is generated by the electronic excitations of glass on a 
time scale comparable to the atomic size divided by the speed of light. Assuming 
instantaneous response, the nonlinear polarization is proportional to a cubic form 
of the field strength. The effect of optical nonlinearities strongly depends on the 




(Al) 
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frequency matching of the field components involved [0] . Consider instead of the full 
electric field E either the positive or negative frequency component oscillating in the 
optical spectral range. To keep the notation simple, we describe these components 
by the symbol E as well. The nonlinear polarization of the medium is only effective 
if P oscillates in a similar spectral range as one of the frequency components of 
E. Consequently, P combines two electric-field components Ei and one complex 
conjugate E^. In isotropic materials such as glass, P is further restricted: isotropy 
implies for the P components [6] 

Pm = ^ {^lEmElEi + K2EiEmEi + K^EiEiE^ (A2) 

where the k denote the material constants of the Kerr nonlinearity. They correspond 
to the x^^^ coefficients in nonlinear optics [6]. In silica, the k constants are identical, 
to a very good approximation, 

Kl = K2 = ^3 = K . (A3) 

In fibers |2], the electric field has two transversal components that we denote by ±. 



Equations (A2) and (A3) give 



P± = 2eo^{ \E±\^E± + '^\E^\^E± + ^E^E*^ ). (A4) 



The first term describes the Self Phase Modulation, the second term the Cross 
Phase Modulation and the third usually is called Four Wave Mixing [2]. The Self 
Phase Modulation does not act across polarizations, whereas the other two terms 
correspond to polarization interactions. For equally polarized fields that carry two 
distinct frequency bands, the Kerr effect acts across these frequencies, because, 
writing 

-ILOat _|_ p ^ — iu>t,t 

we obtain 



E± = Eae-'"^^' + Ei.e-"^"' (A5) 



\E±\''E± = {\Ea\^ + 2\E^\')Eae-''^^' + (l^^l' + 21^,12)^,6-'-"* . (A6) 

The contribution proportional to \Eb\^Ea or \Ea\'^Eb is also called Cross Phase Modu- 
lation [2] . We see that the coupling strength of the Cross Phase Modulation between 
polarizations is 2/3 of the Self Phase Modulation and the Cross Phase Modulation 
within one polarization is twice as strong as the Self Phase Modulation [2]. 

Another cubic nonlinear effect caused by molecular excitations with longer re- 
sponse time, called Stimulated Raman Scattering [2J, contributes to the material 
polarization for pulses below lOOfs duration. Stimulated Raman Scattering is im- 
portant for forming the shape of ultrashort pulses [2], but, it does not act across 
polarizations nor significantly different frequency bands, the cases we are interested 
in; hence we ignore the Raman effect here. 
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A. 2 Waveguides 

Consider the effect of the waveguide on the hght confined in the x and y direction 
and propagating along the fiber in the z direction. We assume that the fiber is 
homogeneous in z and infinitely long, with the Fourier-transformed susceptibility 

X, =X,{^^,x,y) ■ (A7) 

We represent the Fourier-transformed field strengths as 

Eiu,v) = E^iu,z)V^iu,x,y) (A8) 

and require that the fiber modes \J± are eigenfunctions of the transversal part of 
the wave equation for monochromatic light with eigenvalues (3'^{uj), 

V X (V X U±) + (1 + XJ^U±^ = /3|MU± . (A9) 

For single-mode fibers [2j , only one eigenvalue {u) exists for each optical polariza- 
tion ± and frequency u (within a limited frequency range). We normalize the mode 
functions U-t such that the integral of \\J±{u!,x,y)\'^ over the {x,y) plane is unity. 
The optical nonlinearity acts predominantly within a narrow spatial region near the 
fiber core at {x,y) = (0,0). Therefore, we obtain, to a very good approximation, 

[d! + PimE^ = ^ (AlO) 



with the nonlinear polarization (A4). Here we understand k as being averaged over 
the transversal modes (usually causing almost no difference between the two optical 
polarizations ±). 

The eigenvalues {u) of the transversal modes set the effective refractive indices 
n± of the fiber for light pulses E±{t, z) defined by the relation 

/3± = ^a;. (All) 

As a well-known consequence of causality [7] the refractive index must be analytic 
on the upper half plane of complex uj. In the absence of losses within the frequency 
range we are considering, the Fourier-transformed Xg{^) the longitudinal mode 



equation (A9) is real on the real axis and the longitudinal mode equation (A9) is 
Hermitian and positive. Consequently, n± is real on the real axis. Furthermore, 
since the linear susceptibility Xg(t) is real, Xg{^) must be an even function of u, 
which implies that and are even functions of the frequency u. 

When the dielectric structure of the fiber varies over the scale of an optical 
wavelength, the polarization of light becomes an important issue and the refractive 
indices differ for E+{t,z) and E-{t,z), causing birefingence As we have seen, 
microstructured or photonic-crystal fibers |3] allow for some freedom in tailoring the 
dispersion, nonlinearity, and birefringence for specific applications. 
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A. 3 Effective moving medium 

In our case, an intense ultrashort optical pulse interacts with a weak probe field. 
This probe may be caused by an incident continuous wave of light or by the vacuum 
fiuctuations of the electromagnetic field itself [8]. The vacuum fiuctuations are 
carried by modes that behave as weak classical light fields as well. The pulse is 
polarized along one of the eigen-polarizations of the fiber; the probe field may be 
CO- or cross polarized. Due to Stimulated Raman Scattering |2] or the formation 
of optical shocks [2] the pulse will change its shape and velocity; but we assume 
that the intensity profile I{z, t) of the pulse uniformly moves with constant velocity 
u during the interaction with the probe. Since the probe field is weak or in the 
vacuum state we can safely neglect the backaction onto the pulse and the nonlinear 
self interaction of the probe. Since the intensity profile of the pulse is assumed to 
be fixed and given, we focus attention on the probe field. We describe the probe by 
the corresponding component A of the vector potential that generates the electric 
field E and the magnetic field B, with 

E = -dtA , B = d,A. (A12) 

The probe field obeys the wave equation 

{c'dj + c''P^{idt)-dtxdt)A = 0, xoc/(^,t) (A13) 

where x denotes the susceptibility due to the Kerr effect of the pulse on the probe. 



We take the P of Eq. (All) that corresponds to the probe polarization and denote 



the effective refractive index by uq. Equation (A13) shows that the pulse indeed es- 
tablishes an effective moving medium [9] . It is advantageous [2] to use as coordinates 
the retarded time r and the propagation time ( defined as 

T = t--, ( = -, (AM) 
u u 

because in this case the properties of the effective medium depend only on r. Here 
r plays the role of space and C, of time. In this co-moving frame we replace in the 



wave equation (A13) the z and t derivatives by 



and obtain 



dt = dr, d, = -{d^-dr), (A15) 
u 

idc - drfA = dr^n^drA , (A16) 

where the total refractive index n consists of the effective linear index no and the 
contribution due the Kerr effect of the pulse, 

n^ = nl + x. (A17) 

Since x we approximate 

X 

n ^ riQ + 6n , 6n = , (A18) 

2no 

where we can ignore the frequency dependance of uq in x/(2?7,o), which gives Eq. (1) 



of our paper [I]. Note that Eq. (A14) does not describe a Lorentz transformation, 
but the r and ( are of course perfectly valid coordinates; they simply do not belong 
to an inertial system. 
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A. 4 Dispersionless case and metric 



Assume, for simplicity, a dispersionless case where the refractive index no of the 
probe does not depend on the frequency. Note that a horizon inevitably violates this 
condition, because here light comes to a standstill, oscillating at increasingly shorter 
wavelenghts the closer it approaches the horizon. Light waves are dramatically fre- 
quency shifted and thus leave any dispersionless frequency window. However, many 
of the essentials of horizons are still captured within the simplified dispersionless 
model. 



First, we can cast the wave equation (A16) in a relativistic form, introducing a 



relativistic notation [TU] for the coordinates and their derivatives 



X'' 



g- 



(A19) 
(A20) 



and the matrix 

1 -1 
-1 1 — u^n^ /c^ 

that resembles the inverse metric tensor of waves in moving fluids [12]. Adopting 
these definitions and Einstein's summation convention over repeated indices the 



wave equation (A16) appears as 



(A21) 



which is almost the free wave equation in a curved space-time geometry p^. (In 
the case of a constant refractive index the analogy between the moving medium and 
a space-time manifold is perfect [11].) The effective metric tensor g^y is the inverse 
of g^""" [lU]. We obtain 



-1 



(A22) 



In subluminal regions where the velocity c/n of the probe light exceeds the speed of 
the effective medium, i. e. the velocity u of the pulse, the measure of time u^v? j (? — \ 



in the metric (A22) is negative. Here both and dc^ are timelike vectors [T^. In 
superluminal regions, however, c/n is reduced such that u}v? j(? — 1 is positive. A 
horizon, where time stands still, is established where the velocity of light matches 
the speed of the pulse. 



A. 5 Action 

The theory of quantum fields at horizons [131 E] predicts the spontaneous generation 
of particles. In our case, the quantum field is light in dielectric media. The quantum 
theory of light in media at rest has reached a significant level of sophistication [16] , 
because it forms the foundation of quantum optics and, in particular, the quantum 
theory of optical instruments [21 HZ], but quantum light in moving media is much 
less studied [9]. In optical fibers, light is subject to dispersion, which represents 
experimental opportunities on one side, but poses a theoretical challenge on the 



6 



other: we should quantize a field described by a classical wave equation of high order 
in the retarded time. Moreover, strictly speaking, dispersion is always accompanied 
by dissipation, which results in additional quantum fluctuations |16]. Here, however, 
we assume to operate in frequency windows where the absorption is very small. We 
entirely focus on the dispersive properties of the fibre. To deduce the starting point 
of the theory, we begin with the dispersionless case in classical optics and then 
proceed to consider optical dispersion for light quanta. 

The classical wave equation of one-dimensional light propagation in dispersion- 
less media follows from the Principle of Least Action [TO] with the action of the 
electromagnetic field in SI units 



S 



I j ^{n'E^-^B^)dzdt 

j j '-^[n\dtAf-c\d.Af]dzdt 

j j - ^^m^^rA - c\d,Af]udTdC (A23) 



and hence the Lagrangian density 



^ = [AdrU^d^A + c^id^Af] 



(A24) 



In order to include the optical dispersion in the fiber and the effect of the moving 
pulse, we express the refractive index in terms of f3{(jj) and the effective susceptibility 
xij) caused by the pulse, using Eqs. (All ) and ( A17) with ui = idr- We thus propose 
the Lagrangian density 



^ = I [A{c'p\idr) - drxdr)A - c\d,Af] 



(A25) 



In the absence of losses, /3^(ti;) is an even function, as we obtained in Sec. A. 2. We 
write down the Euler- Lagrange equation [TU] for this case 



d{d^A) 



u=0 



^d{d-A) 







(A26) 



and obtain the wave equation (A16). This proves that the Lagrangian density (A25) 
is the correct one. 



A. 6 Quantum field theory 

According to the quantum theory of fields [18] the component A of the vector poten- 
tial is described by an operator A. Since the classical field A is real, the operator A 
must be Hermitian. For finding the dynamics of the quantum field we quantize the 
classical relationship between the field, the canonical momentum density and the 
Hamiltonian: we replace the Poisson bracket between the field A and the momentum 
density jdipc^A) by the fundamental commutator between the quantum field A 
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and the quantized momentum density [18]. We obtain from the Lagrangian (A25) 
the canonical momentum density 



C2 



T: = -eQ—d,A (A27) 
u 

and postulate the equivalent of the standard equal-time commutation relation [181 
[IS] 

[i(C,ri),vr(C,r2)] =-5(n-r2). (A28) 

Til 

We obtain the Hamiltonian 

= -(M)') -^(c'/5'-^rX5r)i)nrfr. (A29) 

One verifies that the Heisenberg equation of the quantum field A is the classical 



wave equation (A16), as we would expect for fields that obey linear field equations. 



A. 7 Mode expansion 

Since the field equation is linear and classical, we represent A as a superposition of 
a complete set of classical modes multiplied by quantum amplitudes d^. The mode 
expansion is Hermitian for a real field such as the electromagnetic field, 

A = J2{AA + Aldl) . (A30) 



The modes Ak obey the classical wave equation (A21) and are subject to the or- 
thonormality relations [131 19] 

{Ak,Aj,,) = 6kk' K,AfeO = (A31) 

with respect to the scalar product 

(Ai, A2) = ^ y iAld,A2 - AAAl) dr . (A32) 
The scalar product is chosen such that it is a conserved quantity for any two solutions 



Ai and A2 of the classical wave equation (A16), 

9^(A,A2) = 0, (A33) 

with a prefactor that turns out to make the commutation relations between the 
mode operators particularly simple and transparent. 

The scalar product serves to identify the quantum amplitudes d^ and a^: the 
amplitude dk belongs to modes Ak with positive norm, whereas the Hermitian con- 
jugate al is the quantum amplitude to modes Al with negative norm, because 

iAl,Al) = -{A,,A2) . (A34) 
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Using the orthonormality relations ( A31 ) we can express the mode operators and 



at as projections of the quantum field A onto the modes A^. and A^ with respect to 



the scalar product (A32), 



{Ak,A) 



■{AlA) 



(A35) 



We obtain from the fundamental commutator ( |A28[ ) and the orthonormality rela- 
tions (A31) of the modes the Bose commutation relations 



[ak,a 



Jkk' 







(A36) 



which justifies the choice of the pref actor in the scalar product (A32). In agree- 
ment with the spin-statistics theorem [TH], light consists of bosons and the quantum 
amplitudes and serve as annihilation and creation operators. 



The expansion (A30) is valid for any orthonormal and complete set of modes. 



Consider stationary modes with frequencies a;^ such that 



dcAk 



-iu'^Ak 



(A37) 



We substitute the mode expansion (A30) in the Hamiltonian (A29) and use the wave 



equation (A16) and the orthonormality relations (A31) to obtain 



H 



a\ak 



1 



(A38) 



Each stationary mode contributes /icj^ to the total energy that also includes the 
vacuum energy. 

The modes with positive norm select the annihilation operators of a quantum 
field, whereas the negative norm modes pick out the creation operators. In other 
words, the norm of the modes determines the particle aspects of the quantum field. 
In the Unruh effect [20], modes with positive norm in the Minkowski space-time 
consist of superpositions of positive and negative norm modes in the frame of an 
accelerated observer [H]. Consequently, the Minkowski vacuum is not the vacuum 
as seen in the accelerated frame. Instead, the accelerated observer perceives the 
Minkowski vacuum as thermal radiation [20] . In the Hawking effect [13], the scat- 
tering of light at the event horizon turns out to mix positive and negative norm 
modes, giving rise to Hawking radiation. 



A. 8 Geometrical optics 

A moving dielectric medium with constant refractive index but nonuniform velocity 
appears to light exactly as an effective space-time geometry [HI [H] . Since a station- 
ary 1 + 1 dimensional geometry is conformally flat [21] a coordinate transformation 
can reduce the wave equation to describing wave propagation in a uniform medium, 
leading to plane- wave solutions [22]. The plane waves appear as phase- modulated 
waves in the original frame. Consequently, in this case, geometrical optics is exact. 
In our case, geometrical optics provides an excellent approximation, because the 
variations of the refractive index are very small. 
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Consider a stationary mode A. We assume that the mode carries a slowly varying 
amplitude A and oscillates with a rapidly changing phase if, 



A = Aexp{i(p) . 



(A39) 



We represent the phase as 



uj{t) dr — ojX 



(A40) 



and obtain from the wave equation (A16) the dispersion relation 



{to - J) 



i\2 



U 



2 2 

n oj 



(A41) 



by neglecting all derivatives of the amplitude A. Here n includes the additional 
susceptibihty x due to the Kerr effect of the pulse according to Eq. 



(A17) 



The dispersion relation has two sets of solutions describing waves that are co- or 
counter-propagating with the pulse in the laboratory frame. Counter-propagating 
waves will experience the pulse as a tiny transient change of the refractive index, 
whereas co-propagating modes may be profoundly affected. Consider the solution 



u 

1 n ) UJ . 

c 



(A42) 



In this case, we obtain outside of the pulse in the laboratory frame cp = n{uj/c)z—u)t, 
which describes light propagating in the positive z direction. Consequently, the 



branch (A42) of the dispersion relation corresponds to co-propagating light waves. 
We also see that u is the frequency of light in the laboratory frame, whereas u' 



is the frequency in the frame co-moving with the pulse. Equation (A42) describes 



how the laboratory-frame and the co-moving frequencies are connected due to the 
Doppler effect. 

In order to find the evolution of the amplitude A, we substitute in the exact scalar 



product (A32) the approximation (A39) with the phase (A40) and the dispersion 
relation (A42). In the limit uj[ — > we obtain 



{A,, A, 



2eoc 
h 



j A?nu exp {iip2 ~ i^Pi)dT , 



(A43) 



which should give S{uj[ — UJ2) according to the normalization (A31). The dominant, 
diverging contribution to this integral, generating the peak of the delta function, 
stems from r ±00 [24]. Hence, for u'l cjg, we replace ip in the integral by if at 
r —>■ ±00 where u does not depend on r anymore. 



2sqc f 

{Ai,A2) = -jr— / A^nu exp [i{uj2 — co'i)r](ir. 



UJ2 - UJi 



duj 

a7' 



^2 - ^'1. 



which gives — uj'2) for 



h 



duj 



(A44) 



(A45) 
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and positive frequencies u in the laboratory frame. Note that positive frequencies 
Lj' in the co-moving frame correspond to negative u in superluminal regions where 
the pulse moves faster than the phase- velocity of the probe light. 

Hamilton's equations [23] determine the trajectories of light rays in the co-moving 
frame, parameterized by the pulse-propagation time (. Here r plays the role of the 



ray's position. Comparing the phase ( A40 ) with the standard structure of the eikonal 



in geometrical optics [4j or the semiclassical wave function in quantum mechanics 
[2l] we notice that —uj plays the role of the conjugate momentum here. Therefore, 
we obtain Hamilton's equations with a different sign than usual [23] . 



r 



duj 



UJ 



(A46) 



We express f in terms of the group index in the laboratory frame. Here the group 
velocity Vg is the derivative of the frequency uj with respect to the wave number 
nuj/c or, equivalently, the inverse of the derivative of nuj/c with respect to uj, which 
gives for the group index c/vg the standard expression [2] 



dn 



ng = n + UJ — 



(A47) 



We obtain from the first of Hamilton's equations (A46) and the Doppler formula 



(A42) 



u 



Ug - 1 



^9 / 

-v' 



— U 



(A48) 



where v'g denotes the difference between the group velocity of the probe Vg and the 
pulse speed u. We see that the velocity f in the co-moving frame (A14) vanishes 



when the Kerr susceptibility x reduces the group velocity c/ug such that it matches 
the speed of the pulse u. Since uj does not vanish here in general, the ray does not 
remain there, but changes direction in the co-moving frame. 

At such a turning point we expect a violation of the validity of geometrical optics 



For example, the amplitude (A45) would diverge here. Geometrical optics is 



an exponentially accurate approximation when 



dT 



dr 



< 1 for T 



27r 



UJ 



(A49) 



as we see from the analogy to the semiclassical approximation in quantum mechanics 
j. Here the cycle T plays the role of the wavelength. We get 



dT _ CjT 
dr ujT 



(A50) 



Consequently, geometrical optics indeed is no longer valid near a turning point where 

(A51) 



c 

u 



This turning point defines a group-velocity horizon where the pulse has slowed down 
the probe such that it matches the speed of the pulse. At this horizon the incident 
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mode is converted into a mode that represents another solution of the dispersion 
relation; a red-or blue-shifted wave, depending on the dispersion and the sign of the 
first derivative of x with respect to r at the group- velocity horizon. White holes 
correspond to increasing x a^nd black holes to decreasing x- White holes blue-shift, 
because incident waves freeze in front of the horizon, oscillating with increasing 
frequency. Black holes red-shift, because they stretch any emerging waves (also 
because black holes are time-reversed white holes). Due to the effective dispersion of 
the fiber, the refractive index changes with frequency. In turn, the dispersion limits 
the frequency shifting by tuning the light out of the grip of the horizon. In particular, 
the dispersion limits the blue-shifting at white-hole horizons to respectable but finite 
frequencies, considering the tiny magnitude of as we discuss in Sec. B2. 

At the event horizons of astrophysical black holes, similar effects are expected 
[25] when, due to the wave-number divergence, the wavelength of light is reduced 
below the Planck length scale ^/hGjc? where G is the gravitational constant. The 
physics beyond the Planck scale is unknown. This trans-Planckian physics should 
regularize the logarithmic phase singularities [13] of modes at the event horizon. 
A numerical study of a simple model of trans-Planckian-type physics [2B] and a 
systematic analytical study [27] indicate that the Hawking effect of the black hole is 
not affected. On the other hand, the quantum radiation of white holes is dominated 
by trans-Planckian physics, because of the extreme blue shift at white-hole horizons. 
It has been predicted [28] that black-hole white-hole pairs could act as black hole 
lasers in a regime of anomalous group-velocity dispersion. From a theoretical point 
of view, trans-Planckian physics regularizes some of the arcane features of quantum 
black holes and gives a more natural picture of the physics behind the Hawking effect 
|27] . In our case, the optical analogue of trans-Planckian physics, optical dispersion, 
is known in principle and turns out to be to the advantage of the experiment. 



A. 9 Classical Hawking effect 

A 'phase-velocity horizon is formed if the pulse has slowed down the probe such that 
its phase velocity is lower than the speed of the pulse. Here an additional effect 
occurs: the spontaneous creation of photon pairs. Hawking radiation. 

In the near ultraviolet around 300nm, the dispersion of microstructured fibers [3] 
is dominated by the bare dispersion of glass where no rapidly grows with frequency 
[2], exceeding the group index cju of the pulse. For such ultraviolet modes, the 



medium moves at superluminal speed. According to the Doppler formula ( A42 ) these 
superluminal modes oscillate with negative frequencies uj' in the co-moving frame for 
positive frequencies lo in the laboratory frame, and vice versa. On the other hand, 
probe modes with phase velocities below c/n oscillate with positive frequencies. 
Therefore, two waves share a given lo' in the co-moving frame, a subluminal wave 
with positive frequency lo in the laboratory frame and a superluminal wave with 



negative cu, see Fig. |A1[ In the case of the astrophysical event horizon, the positive- 
frequency modes correspond to waves outside the horizon that escape into space, 
and the negative-frequency modes to waves beyond the horizon that fall into the 
singularity [14j . The Kerr susceptibility of the pulse may slow down the subluminal 
modes such that the pulse moves at superluminal speed. As we will show, in this 
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Figure Al: Doppler- shifted frequency versus wavelength. The figure shows the co- 
moving frequencies (A42) in units of lO^^Hz for our micro-structured fiber depending on 
the wavelengths A = 27r/tj in nm. To determine the effective refractive index we used the 



dispersion data of our fiber in Eq. (B4), apart from two constants that we fitted to the 



Sellmeier formula for fused silica at short wavelengths [2]- 



case sub- and superluminal modes are partially converted into each other. This 
mode conversion is at the heart of the Hawking effect. 

For simplicity, we consider a single white-hole horizon, not the combination of 
black- and white-hole horizons generated by a moving pulse. We will argue later 
that in practice the white-hole will dominate the Hawking effect, which a-posteriori 
justifies this simplification. Suppose, without loss of generality, that at r = the 
Kerr-reduced phase- velocity of the probe, c/n, matches the group velocity of the 
pulse u. We assume that the mode conversion occurs near this point and expand 
the Kerr susceptibility x as a linear function in r. 



X(^) = Xh + a'r , 



a 



dx 
Or 



(A52) 



The group velocity of the incident probe is much lower than the pulse speed u and 
so both the sub- and the superluminal probe travels from the front of the pulse to 
the back, from negative to positive retarded time r. For a white-hole horizon x 
increases for decreasing retarded time, and so x(0) < 0. 

We proceed similar to Ref. [27] and focus on the conversion region where we 



Fourier-transform with respect to r the wave equation (A16) with the refractive 



index (A17) for stationary waves in the co-moving frame and using the linear ex- 



pansion (A52). The frequency conjugate to r is the laboratory-frame frequency to. 

-id,.,, dt by —iu' and d-r by —iui, denote the Fourier-transformed 



We replace r by —id,^^ 9^ by —iu' and dr by 
vector potential by A, and obtain 



('^o + Xh- ia'dS) ujA 



UJ 
UJ 



2 2 



ujA, 



(A53) 
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This first-order equation has the exact solution 



A 



UJ 



-i(f> 



a" 



UJ 
UJ 



2 2 



duj 



(A54) 



with constant An. We introduce 



u2 a" 



a = —- 



u x(0) 
c 2no 



u dn 
c dr 



(A55) 



in agreement with Eq. (4) of our paper pLj- Note that the phase 
logarithmic contribution, 



contains a 



\nuj + (f)o{uj) 

a' 



(A56) 



This logarithmic asymptotics of the phase will lead to the characteristic mode con- 
version at the group- velocity horizon. In order to see this, we Fourier-transform A 
back to the domain of the retarded time, 



A 



+00 



Ae-'^^^du 



+ 00 



Ao 



UJ 



"duj 



(A57) 



and use the saddle-point approximation, i.e. we quadratically expand the phase 
+ UJT around the stationary points where (9a;(0 + ujt) vanishes and perform the 
integration as Gaussian integrals along the direction of steepest descent. One easily 



verifies that the stationary points are the solutions of the dispersion relation (A41). 



We denote the two solutions by uj± indicating their sign. We obtain for the second 
derivative in the quadratic expansion 



+ u}t) 



n (c 



X(0) UJ \u 



n, 



(A5^ 



The Gaussian integrals at uj± are proportional to the inverse square root of 9^(0 + 



ujt\ We see from Eqs. (A46) and (A48) that they are consistent with the am- 



plitudes (A45) of geometrical optics. Consequently, we obtain a superposition of 



the two waves (A39) that correspond to the two physically- relevant branches of the 
dispersion relation (A41). We denote the positive- frequency wave by A^ and the 
negative-frequency component by A^_. The star indicates that this component re- 
sembles the complex conjugate of a mode, because a mode predominantly contains 
positive laboratory-frame frequencies, according to the normalization (A45). The 
coefficient of A^_ is given by the exponential of the phase integral from the positive 
branch uj^ to the negative frequency uj_ on the complex plane. The amplitude of 
the coefficient is the exponent of the imaginary part of the phase integral, while 
the phase of the coefficient is given by the real part. We can incorporate the phase 
of the superposition coefficient in the prefactor (A45), but not the amplitude. The 



imaginary part of the phase integral comes from the logarithmic term (A56), giving 
-Kuj' joL . Therefore, the relative weight of the negative-frequency component in the 
converted mode is exp(— vrcu'/a'). We thus obtain for r < 



A ~ Z^l'^ [a^ + A*_e-"'^'/"') 



(A59) 
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where Z denotes a constant for given lo' . We determine the physical meaning of Z 
in Sec. AlO, but here we can aheady work out its value by the following procedure: 
consider a wavepacket with co-moving frequencies around lo' that crosses the horizon. 
Suppose that this wavepacket is normalized to unity. After having crossed the 
horizon, the norm of the positive-frequency component is Z, while the negative- 
frequency component has the negative norm — Z exp(— 27ra;'/a'). The sum of the 
two components must give unity, and so 

l-e"2""'/"T\ (A60) 



We represent Z^/^ as cosh^ and obtain from Eq. (A60) that sinh^,^ = Z — \ gives 

Z^l'^ = cosh e , Z^/^g-TTo^'/a' _ gi^j^ ^ _ ^^g^) 

Consequently, the incident wave A± is converted into the superposition A± cosh C, + 
sinh^ when it crosses the horizon from positive to negative r. Hence we obtain 
for this process the mode 

A / ■ > ( 

^ 1 A± cosh^ + At sinhe : r < • ^ ' 



Equation (A62) describes the fate of a classical wave that crosses the horizon. A 
negative-frequency component is generated with weight sinh^ ^ relative to the initial 
wave, but, since cosh^ > 1, the positive-frequency wave has been amplified. The 
mode conversion at the horizon is thus an unusual scattering process where the con- 
certed modes are amplified, at the expense of the energy of the driving mechanism, 
the pulse in our case. (It is also mathematically unusual — the Hawking effect 
corresponds to scattering without turning points in the complex plane.) Wherever 
there is amplification of classical waves, i. e. stimulated emission of waves, there also 
is spontaneous emission of quanta [29] — in the case of horizons. Hawking radiation. 

A. 10 Hawking radiation 

Suppose that no classical probe light is incident; the modes A±ia are in the vacuum 
state. The incident modes are characterized by the asymptotics A± for r > 
while outgoing modes are required to approach A± for r < 0. We perform the 
superposition 

A±o,t = v4±in cosh e - A*^,^ sinh e (A63) 
and see that A±ovLt obeys the asymptotics 

■ f A-t cosh ^ — A* sinh ^ : r>0 /a^,\ 

^--^Ut : r<0 ' ^^6^) 



as required for outgoing modes. The modes (A62) and (A63) describe two sets of 



mode expansions (A30) of one and the same quantum field; for a given u' the sum 



15 



of A-tind-tin and v4^;j^a!j_ijj over the two signs ± of a; must give the corresponding sum 
of v4±outd±out and A*^^^^dl^^^^. Consequently, 

a±in = a±out cosh ^ - a^^^^ sinh ^ (A65) 

and by inversion 

d±out = a±in cosh ^ + a^ij^ sinh ^ . (A66) 

The vacuum state |vac) of the incident field is the eigenstate of the annihilation 
operators a±in with zero eigenvalue (the state that the d-tin annihilate), 

d±in|vac) = 0. (A67) 

To find out whether and how many quanta are spontaneously emitted by the horizon, 
we express the in-coming vacuum in terms of the out-going modes. We denote the 
out-going photon-number eigenstates, the out-going Fock states [9j, by |n+,n_) 
with the integers n±. Using the standard relations for the annihilation and creation 
operators 

a\n) = y/n \n — 1) , a^n) = + 1 |n + 1) , (A68) 
one verifies that d-tin vanishes for the state 



oo 



|vac) = e"'^-'^'/"' n) . (A69) 



n=0 



This is the remarkable result obtained by Hawking [2] for the horizon of the black 
hole. First, it shows that the event horizon spontaneously generates radiation from 
the incident quantum vacuum. Second, the emitted radiation consists of correlated 
photon pairs, each photon on one side is correlated to a partner photon on the other 
side, because they are always produced in pairs. The total quantum state turns 
out to be an Einstein-Podolski- Rosen state [^9], the strongest entangled state for a 
given energy |2Dj. Third, light on either side of the horizon consists of an ensemble 
of photon-number eigenstates with probability z~^e~'^"''^^ . This is a Boltzmann 
distribution of n photons with energies nhuj' and temperature ksT' = ha' /{2ti). 
Consequently, the horizon emits a Planck spectrum of black-body radiation with 
the Hawking temperature (5) of our paper [T]. Fourth, this Planck spectrum is 
consistent with Bekenstein's black-hole thermodynamics [12]: black holes seem to 
have an entropy and a temperature. 

In our case, the spectrum of the emitted quanta is a Planck spectrum for the 
frequencies uo' in the co-moving frame, as long as a phase- velocity horizon exists. We 
performed our analysis for the white-hole horizon, but, since black holes are time- 
reversed white holes, we arrive at the same result for the black hole as well, except 
that the roles of the incident and outgoing modes are reversed. In the laboratory 
frame, the spectrum is given by the dependance of uj' on the laboratory frequency 



u outside of the pulse, i.e. by the dispersion relation (A42) for x = 0. In our 



case, is single-valued for the spectral region where phase- velocity horizons are 



established, see Fig. Al , and so the spectrum of black- and white-hole horizons is 



identical for identical a' . We have shown in Eq. (8) that this mapping from uj to 
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uj' amounts to a re-definition of tlie Hawking temperature: in tfie laboratory frame 
ksT is given by tlie logaritfimic derivative of x the horizon; it is independent 
of the magnitude of the Kerr susceptibihty, as long as a phase-velocity horizon is 
established. The particle-production rate depends only on the sharpness of the pulse. 
This important feature makes the experimental observation of Hawking radiation in 
optical fibers feasible using few-cycle pulses [HI] . 



A. 11 Optical shock 

Another feature increases the Hawking temperature further: the formation of an 
optical shock-front at the trailing end of the pulse. For a sufficiently intense pulse, 
the Kerr effect strongly influences its shape; it counteracts the dispersion of the 
pulse, forming a soliton [2] and it may lead to self-steepening [32]. Consider the 
propagation of the pulse itself. We describe its electric field E by the envelope £ 
and the phase at the carrier frequency ojq, 

E = S exp (i(3oz — iujQt) , /?o = ''^0(1^0 )— , -^=1^1^- (A70) 

c 

Assuming that the envelope varies over longer scales than an optical cycle, we ap- 



tion (A4) by 



proximate in the wave equation (AlO) the differential kernel and the Kerr polariza- 



dlP ^ -2ieoKUoE {-ituol + WJ) . (A71) 

The dominant nonlinear term in the resulting wave equation is proportional to the 
intensity /; the term proportional to dtl only becomes important for sharp features 
and, as we will see, leads to the formation of an optical shock. 

First we ignore the shock term and approximate 13{(jj) around the carrier fre- 
quency ojq by a quadratic polynomial, 

(3^(3o + (3i{uj-uj^) + ^{uj- ujof ■ (A72) 

Here /3i describes the inverse of the group velocity u and (32 the group-velocity 
dispersion. Substituting idt for uj ~ujq wherever it acts on the envelope £ we obtain 
the usual nonlinear Schrodinger equation of Nonlinear Fiber Optics [2] 

+ ^d) £ - ^d^£ + -^\£\^£ = 0. (A73) 

This equation is integrable by the Inverse Scattering Method [30] . For anomalous 
group-velocity dispersion, where P2 < 0, it has the fundamental soliton solution 

in terms of the retarded time r and the propagation time C, according to the definition 



^ = £:osech( -Jexp^z^J , ^0 = ^^ (A74) 



(A14). The constant Tq describes the duration of the soliton. 
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Consider now shock formation. For simplicity, we ignore the group-velocity dis- 
persion and arrive at the nonlinear wave equation 



d^S + j£{~zujoI + 3drI) = 0, 1 = ^^. (A75) 

C PoPl 

We represent £ in terms of the amplitude vT and a phase, and obtain that the 
evolution of the phase is completely determined by the amplitude, whereas the 
intensity obeys the transport equation 

{dt; + Q-fldr) 1 = 0. (A76) 

The intensity is transported with velocity 67/ where the retarded time r plays the 
role of space and the propagation time ( the role of time. In a (r, () space-time 
diagram, the initial intensity profile /o(to) is thus transported along a line with 
steepness G'jI^tq). Hence we write down the solution as 

/ = /o(ro), r-67/o(ro)C = T-o. (A77) 

The relationship between r, ( and the initial intensity profile implicitly determines 
tq. At some r and (, two lines of transported intensity, one belonging to Tq and the 
other to Tq, may cross, developing a discontinuity, a shock. Here 

C ~ To "To 

The shock time ( is smallest for the largest negative derivative of the initial intensity 
profile. So the steepest point at the trailing end of the initial pulse is the first to 
form a shock. Although we ignored the dispersion and Stimulated Raman Scattering 
[2] in this simple theory of optical shocks [32] ultrashort pulses still form sharp 
features at their trailing end, features that may easily become comparable to the 
carrier wavelength, as we and others have seen in numerical simulations using the 
method of Ref. [31]. Therefore, the white-hole horizon will dominate the Hawking 
radiation of an optical pulse and the radiation is likely to become strong enough to 
be detectable. 
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B Experiment 



In this appendix we describe the experimental observation of frequency shifting of 
hght at the group velocity horizon. To our knowledge, such an experiment has 
never been carried out before; the cases closest to our scheme are demonstrations 
of pulse trapping [321 ES] where the dynamics are dominated by the Raman effect 
or pulse compression in a fiber grating (optical push broom) j|37j. Based on the 
theory of Appendix A, we also derive mathematical expressions for the amount of 
blue shifting, for the spectral shape, and for estimating the efficiency of this process. 
We discuss the experimental proceedings and findings and compare them with the 
theory. 



B.l Dispersion 

The creation of artificial event horizons in optical fibers critically depends on the op- 
tical properties of these fibers. These properties are described here and summarized 



in Table B.l To create an artificial event horizon in our scheme, an intense optical 
pulse has to be formed inside the fiber. Optical solitons [3HI Ej offer a unique pos- 
sibility for nondispersive stable pulses in fibers. These can be ultrashort, allowing 
for very high peak powers to drive the nonlinearity of the fiber. Bright solitons only 
exist for anomalous group velocity dispersion [2]. Microstructured fibers |3] have an 
arrangement of holes close to the fiber core along the fiber. In the simplest picture, 
the holes lower the local refractive index in the transverse plane of the fiber, leading 
to substantially larger index variations compared to conventional fibers. Hence a 
very wide range of transverse refractive index profiles can be engineered. There are 
various designs for the shape and location of the holes, leading to a range of different 
effective dispersions and giving rise to a variety of applications [S]. In particular, the 
anomalous dispersion required for solitons can be generated at wavelengths reaching 
the visible. 

The dispersion parameter D of optical fibers is defined as the change of group de- 
lay per wavelength change and propagation length. Its units are usually ps/(nmkm). 
Since the group delay per propagation length is given by n^/c and ng/ c=d(3 / doj, we 
have [2] 

oXouj uj 

The group- velocity dispersion is often also characterized by the second derivative of 
(5 with respect to the frequency tu. 

The group- velocity dispersion is normal for positive j32 and negative D, and anoma- 
lous for negative (32 and positive D. 

For the creation of a horizon we chose a commercial microstructured fiber, model 



NL-PM-750B by Crystal Fiber A/S. Figure Bl shows the dispersion of the particular 
fiber sample we used. The red curve is a measurement provided by Crystal Fiber; the 
dotted line was measured for our fiber sample by Alexander Podlipensky and Philip 
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Figure Bl: Two measurements of the dispersion parameter D for the fiber used in the 
experiments. Red hne: company data; dotted Hne: result of Alexander Podlipensky and 
Philip Russell at the Max Planck Research Group in Optics, Information and Photonics 
in Erlangen, Germany. If the total shaded area vanishes, the two wavelengths at either 
end are group velocity matched. 



Russell at the Max Planck Research Group in Optics, Information and Photonics 
in Erlangen, Germany. The fiber dispersion is anomalous between 740 nm and 
1235 nm wavelength and normal otherwise. Thus solitons can be created using 
ultrashort pulses from modelocked Ti:Sapphire lasers. Light that would probe the 
horizon and experience blue shifting as a result, will have to be slowed down by the 
Kerr effect of the pulse such that its group velocity matches the speed u of the pulse. 
The Kerr susceptibility is small (we give an estimate in Sec. B2), and so the initial 
group velocity of the probe should be only slightly higher than u. Integrating Eq. 



(Bl ) we obtain 

A 



Dd\ = /?i(A) - /?i(Ao) = — ^ - — ^ , vg{\^) = u . (B3) 
Ao ""gl^J Vg[^Ao) 

Here Aq and A denote the center wavelengths of the pulse and the probe light, 
respectively. Therefore, the probe light travels at the speed of the pulse if the integral 
of D vanishes, as illustrated by the shaded areas in Fig. |Bl[ This probe wavelength is 
called the group-velocity-matched wavelength Am and the corresponding frequency 
iOm the group- velocity-matched frequency. For a pulse carrier-wavelength of 800 nm 
and the fiber used here we obtain Am ~ 1500 nm. This value of Am is useful, because 
on the one hand it is a standard wavelength for lasers and optical equipment made 
for fiber-optical telecommunications and on the other hand it is clearly separated 
from our broadband pulsed light. 

The dispersion D essentially describes the effective group index of the fiber as 



a function of wavelength (or frequency). Integrating Eq. (Bl) twice we obtain the 
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Property 


Crystal riber \ brlangen 


Dispersion Dq 


28\36ps/(nm km) 


Dispersion 


-150 \ - 180ps/(nmkm) 


Third order dispersion dD„Jd\ 


-0.6 \ - 0.75ps/(nm2km) 


Dispersion /32(Ao) 


-9.5 \ - 12psVkm 


Dispersion /32(Am) 


180\210psVkm 


Group velocity-matched wavelength Am 


1508\1494 nm 


Birefringence Ar?,o\Ar?,m 


7.5\5.7 X 10-^ 


Nonlinearity 7 (780 nm) 


0.1 W-^ii-i 


Fiber length L 


1.5 m 



Table 1: Properties of fiber NL-PM-750B. Dispersion data according to Crystal Fiber 
\ Alexander Podlipensky and Philip Russell, Max Planck Research Group in Optics, In- 
formation and Photonics in Erlangen. Nonlinearity according to Crystal Fiber. The 
birefringence An and the fiber length L were measured by the authors. The symbols used 
are defined in the text. 



effective wavenumber (All) 



(5 = /3(cuo) + ^ ~ ^0 _^ / DdXduj = n- 

= Piuo) + "^^^^ + [ [ /32dudu, (B4) 



u 



where n is the linear effective refractive index of the fiber and Uq and u denote the 
carrier frequency of pulse and probe, respectively. So, in addition to the measured 
dispersion curve, only two constants determine (3: the pulse group velocity u at ujq 
and P{lvo). However, in what follows the most relevant parameters for describing 
our experiment are independent of these constants. 

In general, the two polarization modes of the fiber have slightly different prop- 
agation constants /?. This birefringence creates a refractive index difference of An 
between the polarization modes. In the fiber chosen in this experiment, the holes 
form a hexagonal pattern that is slightly distorted to break the hexagonal symme- 
try. Our fiber exhibits strong birefringence An of a few times 10""^. This leads to 
non-negligible changes in the group velocity. 



B.2 Frequency shifts 

Consider the frequency shifts at a group velocity horizon. During the pulse-probe 
interaction, the co-moving frequency u' is a conserved quantity and so the probe 
frequency u follows a contour line of uj' as a function of the nonlinear susceptibility 



X induced by the pulse, see Fig. |B2[ The maximal x experienced by the probe is 
proportional to the maximal nonlinear susceptibility xo experienced by the pulse: 
assuming perfect mode overlap of pulse and probe, Xmax reaches 2xo when the probe 
and the pulse are co-polarized and 2xo/3 when they are cross-polarized, see Sec. Al. 



21 



If the pulse is a soliton, we obtain from Eq. (A74) the peak susceptibihty 



Xo 



2nocXoDo 2noc|/52(Ao)| 



(B5) 



where Dq denotes the dispersion parameter at the carrier wavelength Aq. For exam- 



ple, for a soliton (A74) at Ao=800nm whose full width at half maximum (FWHM) 
is 70fs (corresponding to To=40fs), for no=1.5, Dm=30ps/(nmkm) the peak sus- 
ceptibility Xo is as low as 2 X 10^^. Nevertheless, we show that this small variation 
in the optical properties is sufficient to generate a significant wavelength shift at the 
horizon. 
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Figure B2: Doppler contours. The pulse shifts the laboratory frequency uj (or the 



wavelength A) along the contour lines of uj' as the function (B6) of the refractive-index 
change 6n = x/(2no). For a sufficiently intense pulse 5n reaches the top of a contour. In 
this case the probe light completes an arch on the diagram while leaving the pulse; it is 
red- or blue-shifted, depending on its initial frequency. 



We obtain the contours of lu' from the Doppler formula (A42). We use relations 



(|A18|) and (|B4|), but integrate from the group velocity-matching point, 

X u 



UJ 



uj'^ - u I I (32 duo duo 

I UJm ^ 

u 

I TTU 
Am 



■ UJ 



Am, 



2^0 c 
X u 
2nQ c 

_ X_ 
no 



■UJr. 



(B6) 
(B7) 



The contours of uo' do not depend on uj'^ nor on the scaling factor T^u/Xm- Because 
Dm = -D(Am) < 0, they form inverted parabolas with a maximum at A^, for the 
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corresponding Xmax- They intersect the axis of zero x the incident and the 
emerging wavelengths. Here \Dm\c{\ — Xm)"^ l^m equals Xmax/'^O) and so we get 



A = Am ± (5A , 5\ 



AmXn 



\Dm\ noc 



(B8) 



Using again that the pulse is a soliton, we obtain 



SX 



2r Ar 



Tq^UJqUJ 

n 



Dn 



V2rAmAo 

To UJm 



/?2(Ao) 



(B9) 



with r=2 for co-polarized and r=2/3 for cross-polarized pulse and probe light. Ac- 
cording to Fig. |B2| the probe light can maximally be wavelength-shifted from +6X to 
—6X over the range 26 X. For the soliton mentioned above the group velocity disper- 
sion Do is about 30ps/(nmkm). Using Am ~ 1500 nm and ^ — 160 ps/(nmkm), 
the wavelength shift 26X is 20 nm in the co-polarized case and 25A=12nm in the 
cross-polarized case. 

We also derive a simple estimate of the efficiency of the frequency shifting from 
the dispersion data. The probe light that is colliding with the pulse undergoes 
frequency conversion at the horizon. However, because the group velocities of the 
probe Vg and of the pulse u are similar, only a small fraction of the total probe light 
can be converted within the finite length of the fiber. The pulse and the slightly 
faster probe hght travel through the fiber in t=L/u and tp=L/vg with t > tp. The 
time difference multiplied with the probe power Pprobe is the energy Ecoii converted 
by pulse collision: Ecoii=PprohcL {1/u — 1/vg). Therefore, the fraction t] of probe 
power that is frequency converted is 



1] = UrcpL{l/U - l/Vg) 



L\Dm\SX, 



(BIO) 



where u^ep is the repetition rate of the pulses and l/vg=dp/duj ^ 1/u + Dm6X was 
used. For L=1.5m and z/rep=80MHz the maximal conversion efficiency r] is on the 
order of 10~^. 



B.3 Experimental results 



The experiment is arranged as displayed in Fig. B3 A modelocked Ti:Sapphire 
laser (Mai Tai, Spectra Physics) delivers 70-fs pulses (FWHM) in the near infrared 
(NIR). These linearly polarized pulses are coupled to either one of the principal 
axes of the microstructured fiber of length L=1.5 m. The polarization is rotated by 
a half-wave plate. Note that the polarizing beam splitter (PBS) only acts on the 
probe light. At the fiber output temporal autocorrelation traces and spectra are 
taken to determine the pulse energy necessary to create a fundamental soliton. For 
the center wavelength of 803 nm, a dispersion Do=30 ps/ (nmkm) and a nonlinearity 
7 of 0.1 W~^m~^, 70-fs solitons are generated at 5pJ pulse energy corresponding to 
400 /iW average power for the repetition rate i'^cp=80 MHz. 

The output pulse length equalled the 70-fs input pulse length at an input power 
of approximately 320 /^W. This indicates that a soliton has formed. The observed 
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cw probe laser (IR) 



pulsed Ti:Sapphire 
laser (NIR) 




Figure B3: Light from a continuous- wave infrared (IR) laser is reflected off a diffraction 
grating (G) to suppress fluorescence. The beam is steered into a microstructured fiber 
(MF) by a polarizing beam splitter (PBS), a half wave plate (HWP), and a coupling 
lens (L). Near-infrared (NIR) pulses are launched as well. After the fiber the pulses are 
removed by a dichroic filter (DF) and the probe spectrum is taken through a standard 
single mode fiber (SMF). 



power of 320 fiW in comparison with the predicted power of 400 /iW illustrates the 
uncertainty in the actual fiber dispersion and nonlinearity. The observed Raman- 
induced soliton self- frequency shift [21 139] was < 4 nm. Note that this shift deceler- 
ates the pulse and hence is changing the group velocity-matched wavelength Am in 



the infrared (IR). To calculate how much is shifted, we use Eq. (B3), replacing 



Ao and A with Aq + 5Ao and A^ + ^A^ and linearize. In this way we get 

SXr = ^SXo. (Bll) 



For the dispersion data shown in Fig. Bl, a wavelength change of 4nm of the pulse 



changes A^ by ^A^ = — 0.75nm. Since the probe light is wavelength-shifted sym- 
metrically around A^, there is a change of the wavelength shift of up to — 1.5nm. 

The probe light is derived from a tunable external grating diode laser (Lynx 
Series, Sacher Lasertechnik) . It delivers up to 20 mW of continuous-wave light, 
tunable from 1460 to 1540 nm. The probe light is reflected off a diffraction grating 
to reduce fluorescence emitted near lasing bandwidth. With another half-wave plate 
the probe light is coupled into the fiber onto one of the principal axes. Depending 
on wavelength, 100 to 600 /xW of probe power were coupled through the fiber. After 
the fiber we use a dichroic optic to filter out all of the pulse light and couple the IR 
light into a single-mode fiber connected to an optical spectrum analyzer. 



Figure B4 shows a typical output spectrum. This spectrum was taken with pulse 



and probe aligned to the slow axis of the fiber. At A=1506nm the diode-laser input 
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Figure B4: Spectrum of the blue-shifted hght for initial probe of wavelength 1506nm. 
Pulse and probe are co-polarized along the slow axis of our fiber. Traces with (green) and 
without (black) probe light are shown on a logarithmic scale. They are subtracted on a 
linear scale to obtain the normalized signal (red) displayed in parts per million. 



line is visible as a strong signal. From A=1502nm to A=1510nm we detect residual 
weak spontaneous emission from the laser that was not completely eliminated by 
the diffraction grating. Traces with and without pulses present in the fiber are taken 
and subtracted, leading to the signal displayed on a linear scale (red color). The 
signal is normalized by the amount of probe power and by the resolution bandwidth 
of 0.5 nm and is given in parts per million (ppm). With the pulses present, a clear 
peak appears on the blue side of the input probe light near 1493 nm. Since the blue- 
shifted light is generated from the part of the probe light that overlapped with the 
pulse during fiber propagation, it constitutes itself a pulse of finite length. Hence, 
this length is determined by the relative group velocity of probe light and the pulse, 
see for example Eq. (BIO). In turn, the unshifted probe light is partially depleted, 
forming a gap in intensity. These features lead to a spectral broadening of both the 
shifted and unshifted probe light by a few nanometers. 

From the measurements shown in Fig. B4, the efficiency of the blue-shifting is 
1.1 X 10~^, less than the estimated 10~^. This indicates that a significant part 
of the probe light tunnels through the pulse; the pulse is too short to establish a 
nearly perfect barrier. In the tunneling region of the pulse the laboratory frequency 
uo is purely imaginary. In order to estimate the maximal imaginary part of uo we 
consider the extreme case where the initial frequency of the probe reaches the group- 
velocity-matched frequency uj^ characterized by uj' = u'^. We solve Eq. (B6) for u 
and obtain 



lmu=J^^= , . (B12) 



Assuming x ^ 2 x 10 ^ at the soliton peak, no ~ 1.5, Am = 1500 nm and Dr, 
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— 160 ps/(nmkm) the imaginary part of uj reaches about 5THz. This is insufficient 
to significantly suppress tunnehng through a 70-fs pulse, because the product of 
Ima; and Tq is much smaller than unity. For longer pulses we would expect perfect 
frequency conversion at the horizon. 

Increasing the probe wavelength further away from is shifting light further to 
the blue side of the spectrum, because the wavelength shifts symmetrically around 



the group velocity-matched wavelength, according to Eq. (B8) and Fig. B2 Figure 



B5 displays the spectra of shifted light for three detunings of the probe light from 



the group velocity-matched wavelength Am- As expected, the spectra move towards 
shorter wavelengths by the same amount as the probe laser was tuned towards longer 
wavelengths. 




Figure B5: Spectra for different input probe wavelengths. Since the probe mode is 
mirrored around the group- velocity-matched wavelength Am, increasing probe wavelengths 



experience increasing blue shifting, as is also illustrated by the contours of Fig. B2 



Figure B6 shows how the signal strength, the spectrum integrated over the signal 
peak, evolves with increasing probe power. A clear linear dependence is evident in 
agreement with our theoretical model and the superposition principle for the probe 
light. The figure illustrates that the probe indeed is a probe, not infiuencing the 
pulses via nonlinear effects. 

Changing the input polarizations changes the group velocities of pulse and probe 
and therefore the group velocity-matched wavelength A^, shifts by an amount 5Xm- 
If we change for example the pump polarization from the fast to the slow axis, the 
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Figure B6: Integrated blue-shifted spectrum versus incident probe power. The solid line 
is a fit assuming proportionality. 



inverse group velocity (3i{Xo) increases as rigo/c is replaced by [ugQ + Ano)/c. To 
maintain group velocity matching, /3i(Am) has to change accordingly by Auq/c. We 



use Eq. (B3), linearizing around Am, and get 

Ano/c ^ DmSXm . (B13) 

For Z)m=~160ps/(nmkm) and Ano=7.5 x 10^^ we obtain 6Xm ~ — 16nm ^ 26X. 
This means that, when changing polarizations, the probe laser has to be retuned to 
a wavelength were frequency shifting can be observed. 



Figure |B7| shows spectra for all four different polarization combinations. As 
expected, the group velocity-matched wavelength changes. Note that there also is 
a difference in for the two co-polarized cases, indicating small changes in the 
dispersion profile for the two polarization dispersion of the birefringence. 

B.4 Calculation of spectra 

Here we derive a functional expression for the spectra of the frequency-shifted probe 
light, applying the theory developed in Appendix A. We focus on the part of the 
incident light that interacts with the pulse and ignore the component that tunnels 
through, because the tunneled light does not contribute to the observed spectrum. 



We represent the relevant probe light A as a superposition (A30) of stationary 
modes Aji and on the right or left side of the horizon that are characterized by 
the frequency a;' in the co-moving frame. We replace the mode operators by classical 
amplitudes aji and and focus on the positive-frequency component of A. In this 
way we obtain 

A = j (anAn + glAl) duo' . (B14) 

For the mode functions Aji and A^ we use the theory of geometrical optics in 
moving media developed in Sec. A. 8. We assume that the pulse is infinitely short 
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Figure B7: Blue-shifted spectra for all four polarization combinations. Co-polarized 
spectra on the slow and fast axis in (a) and (b) and cross-polarized spectra with pulses 
slow (c) and fast (d). Group velocity-matched wavelengths are (a) 1499.5 nm (b) 1503.2 nm 
(c) 1486.4 nm and (d) 1513.3 nm. 



in comparison with variations of the probe light; and write down the modes as 

Ar = e(r) (Ae~'"^" + M2e-*"^")e-*"'S 

Al = 0(-r) (Ae"*^^" -Mse-*'^^") e"*'^'^. (B15) 

Here Ui denotes the laboratory frequency with respect to u' and uj2 the blue-shifted 
laboratory frequency. The factors of ±i describe the 11/2 phase shifts at turning 
points [21]. Their sign depends on the sign of the frequency change and the side of 
the horizon. We assume that ui and U2 are sufficiently close to the group velocity- 
matched frequency Um such that we can use the quadratic approximation ( B6 ) for 
uj' outside the pulse where x vanishes, 

^' = ^'m - 2/^2 - ^mf , (B16) 

which implies that the two solutions lox and 102 for a given lo' are symmetric around 
the group velocity-matched frequency cUm, such that 

<^2 = + {^m - ^\) = '2idm - ^1 ■ (BIT) 
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The amplitudes Ai and A2 are given by the expression (A45). Since nu does not 
vary much with frequency and duji/ duj' =—duj2/ duj\ we approximate 



Ai^ A2 = Aq 



duji 

a7 



(BI8) 



Finally, we write the mode expansion (B14) as an infinite integral over uji, assuming 



that contributions outside the physically relevant range of uoi are negligible, 



A 



+00 



[an{ui)AR + aiMAi) dui 



(B19) 



Formulas (B15 B19) specify our theoretical model. 



In our experiment, we measure the modulus squared of the Fourier transform of 
A with respect to the laboratory time t at the end of the fiber z=L. This is identical 
to the modulus squared of the Fourier transform A with respect to the retarded time 
T=t — z/u dX C,=L/u. Using the standard relations 







+ 'k5{uj) 



(B20) 



we obtain the result 



A 



271 



^0 



A(r)e^'^" dr 



aniuj) + aL{uj) 



-iu)'C 



+ iAi 



aii{2ujr, 



UJ 



UJ 



+00 



27r(u; — iUi) 



271(00 — UJ2] 



duJi (B21) 



where the integral is understood as a Principal Value Integral [IQ] through the poles 



where Ui or UJ2 go through u. In the first line of Eq. (B21 ) cj' is a function of u and 



in the integral in the second line u' is understood to be a function of Ui. 

At the entrance of the fiber, the initial probe Ap is a continuous wave with 
frequency Up and an amplitude we denote as 2Ao. The coefficients qr and ql that 
describe this situation are given by the expressions 



ttR = S{uJi - Up) + 



7C{UJI - Up) 



ai = S{uji - Up) - 



7t{ui - Up) 



(B22) 



as one verifies by the following procedure: we substitute the coefficients in the 



integral (B19), extract the contribution of the delta function and apply Cauchy's 
lO] for the remaining integral. For this, we close the integration 



Residue Theorem 

contours on the complex half planes where the integrand exponentially decreases. 
For the modes on the right-hand side of the horizon, where r > 0, we chose the lower 
half plane for the exp(— zcjit) term and the upper half plane for exp(— iu;2''"), in view 
of the relationship ( B17[ ); on the left-hand side we take the opposite planes. Since 
we integrate through the poles we obtain half of the residue, similar to the derivation 
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of Hilbert transformations [10] or Kramers- Kronig relations [7]. The result of this 



calculation is the incident plane wave 2Aoexp{ 



which justifies the mode 



coefficients (B22). 



In order to calculate the spectrum, we substitute the mode coefficients (B22) 



into Eq. (B21). The delta functions immediately generate contributions to the 



spectrum; we focus on the calculation of the remaining Principal Value Integral 

with Gaussian exp(— icu'i^) given by Eq. (B16). We 

-1 



through [uj 

represent {uj — ujp){uj — uji)^^{uji — ujp 
Hilbert transform of a Gaussian |i41j . 



as [ui — Up 



[Ui - UJ 



and use the 



+00 g-e 



tt{x - 



e ^ erfi(x) , erfi(x) 



^d^. (B23) 



In this way we find that the spectral field A consists of the sum of the delta peak 
Ao6{lj — Up) exp(— zcijpC) and the contribution that describes the frequency shifting 
of the probe, 



7!-{lJ + LJp- 2LJm 



exp {iq{uj - Umf) - exp {iq{ujp - UmY) 



As 



TC{UJ — UJp) 

where we ignored the unimportant overall phase of uj'^C, and used the abbreviation 



G ( (t^ - t^m) ] - G [ {uJp - UJm) ) (B24) 



TTc\Dm\L 
ujI 



(B25) 



So, up to an overall phase, the spectral field A is given in terms of experimentally 
accessible parameters, the fiber length L, the dispersion and the group velocity- 
matched frequency uJm- 

As expected, the spectrum gets narrower around the blue-shifted and probe fre- 
quencies 2ujm — ^jJp and uip with increasing propagation distance in the fiber, because 
both the converted and depleted components of the probe light form longer pulses 
for longer interaction times. However, the carrier frequency of the pulse is grad- 
ually red-shifted due to the soliton self-frequency shift (SFS) [22] • According to 



Eq. (Bll), this leads to a shift 5\r in the group velocity-matched wavelength 
of about — 0.75nm along the fiber. The blue-shifted light is created at decreasing 
wavelengths and also the part of the pulse that interacts with the probe adiabat- 
ically follows. In the signal spectrum, both the group-velocity-matched frequency 
UJm and the probe frequency ujp appear to be shifted by 5uJr = —{6Xr/ Xm)^m and 
is replaced by ujg. Since initially the spectrum is broad, we incorporate the SFS 
effect in a phenomenological form in our formula by a reduced efficiency rir for the 
blue-shifted part of the spectrum as 



As 



VrAo 



Tt{uJ + UJs- 2uJm) 



exp [iq{uj - uJm)'^) - exp {iq{uJs - Um)^) 



n{uj - UJs) 



G [y^^{uj - UJm)) - G [y^^{uJs - UJm)) . (B26) 
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Figure |B8| shows the fit of the observed spectrum of Fig. |B4| with the theoretical 



curve (B26). As fitting parameters we used the overall amplitude Aq, the shifted 
probe frequency Us and group-velocity-matched frequency Um (in terms of the cor- 
responding wavelengths), the dispersion Dm and r]r. We obtain a very good fit for 
^0 = 4.1 X lO^^VWs, Xs = 1505.31 nm, = 1499.38 nm, = -187ps/(nmkm) 
and ?7r = 0.80. The shift in Xp is consistent with the effect (Bll) of the soliton 

and Xr 



self-frequency shift. The fitted values for Dm and agree with the independently 
measured dispersion and the group-velocity-matched frequency calculated from the 



dispersion curve of Fig. Bl 



In conclusion, we have shown that light was blue-shifted by a near group velocity- 
matched pulse. The measured data was explained by the presence of an optical group 
velocity horizon inside the fiber. A very good agreement between theory and ex- 
periment was achieved. The blue shifting corresponds to the optical analogue of 
trans-Planckian frequency shifts in astrophysics [25j. In this way, we have demon- 
strated classical optical effects of the event horizon in our analogue system, a first 
step towards tabletop astrophysics 
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Figure B8: Fit with theory. The red curve shows the normalized difference data of Fig. 



B4 and the thin black line the fit with the theoretical expression (B26|. 
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